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I.  Experimental  Dynamic  Response  of  a 
Viscoplastic  Beam  Segment 

R,  DiMarcello  and  V.  H.  Neubert 


Background 

The  writers  and  associates  have  conducted  a  series  of  tests  of 
viscoplastic  beams  under  dynamic  loading.  The  relation  of  this  work  to 
that  of  others  in  the  literature  has  been  summarized  in  past  reports ^ 
and  will  not  be  reported  here.  However,  enough  of  the  local  work  will 

be  reviewed  to  show  the  motive  for  the  study  reported  herein. 

(2) 

Stanovsky  carried  out  experiments  on  cantilever  steel  beams  of 
rectangular  cross-section.  The  beams  earned  a  tip  mass  and  were 
loaded  by  a  mechanical  impactor  which  struck  the  tip  mass.  A  pressure 
transducer  between  the  mass  and  the  impactor  measured  the  applied 
force,  which  was  used  to  predict  the  measured  strain,  One  of  the 
primary  conclusions  of  Stanovsky's  work  was  chat,  for  mild  steel,  the 
experimental  strain  lagged  behind  the  strain  predicted  using  a  bi¬ 
linear  stress-strain  relationship,  This  led  to  the  study  of  the 
following  moment-curvature  relationship 

EIk»M  +  C(M-Mj  (1) 

st 

(3) 

by  Brown  .  Equation  (1)  is  similar  in  form  to  a  stress-strain 
suggested  by  Malvern ^  and  studied  by  Plass^  and  others-  A 
derivation  cf  Equation  (1)  starting  with  the  stress-strain  relationship 
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is  given  in  Part  II  of  this  report.  Brown  performed  experiments  on 
cantilever  beams  in  a  shock  machine.  Moments  were  calculated  based  on 
measured  accelerations  of  the  tip  masses-  Curvature  was  deduced  from 
strain  measurements-  Brown  suggested  values  of  C  for  Equation  (1)  in 
the  range  of  C  =  500  to  3000  for  a  1"  x  1/8"  beam.  Subsequently  tests 
involving  three  parallel  beams  carrying  a  common  tip  mass  were  carried 
out  on  the  shock  machine  at  the  Naval  Research  Laboratory.  Again  the 
moments  were  computed  from  measured  accelerations  of  the  tip  mass. 
Static  tests  of  beams  of  the  same  material  were  made*  Vogel^*^  pre¬ 
dicted  the  dynamic  response  from  measured  base  acceleration  and  found 
that  Equation  (1)  a  satisfactory  moment-curvature  relationship  using 
values  of  C  in  approximately  the  same  range  as  proposed  by  Brown-  One 
unexpected  result  of  these  tests  at  the  Naval  Research  Laboratory  was 
that  the  measured  dynamic  moment  was  less  than  the  measured  static 
moment  just  after  yielding  of  the  outer  beam  fibers-  It  was  not 
determined  whether  this  was  a  property  of  che  material  or  due  to  errors 
in  measurements-  The  main  purpose  of  the  present  experiments  was  to 
study  the  same  beam  cross-section,  but  to  use  a  different  experimental 
method  for  deducing  the  applied  moment 

Description  of  the  Apparacus 

The  apparatus  is  shown  m  Figure  )  The  beam  wa s  .1  inch  wide, 

1/“  inch  deep,  and  18  inches  long.  It  carried  a  steel  tip  mass  bolted 
on  at  each  end  There  were  twc  supports  which  weie  designed  to  act 
approximately  as  simple  hinge-*-ype  supp.rts.  The  experimental  segment 
was  a  machined-dcwn  legion  at  the  -sn'et-une  1'8  inch  thick,  1  inch 
wide,  and  L  inch  long.  The  apparatus  had  two  features  that  were 
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significantly  different  from  those  of  previous  tests:  (1)  the  test 
segment  was  in  a  region  to  which  the  shear  should  have  been  a  minimum 
and  (2)  the  moment  was  deduced  from  strain  gauges  at  locations  1  and  4 
where  the  beam  material  remained  essentially  elastic. 

The  beam  supports  were  bolted  to  a  base  which  was  in  turn  bolted 
tc  the  table  of  the  IMPAC  shock  machine,  the  same  shock  machine  used  by 
Brown  The  acceleration  pulse  provided  at  the  shock  machine  table  is 
approximately  a  half-sine  wave  with  a  20  millisecond  duration. 

Strains  were  measured  using  Micromeasurement  foil  gauges  1/8  inch 
long  Ellis  BAM-1  bridge  amplifiers  were  used  with  one  arm  active  and 
the  signals  were  photographed  cn  a  Tektronix  oscilloscope  (Figure  2). 

Stitu.  tests  were  performed  using  the  same  beam  and  support 
arrangement  cn  the  table  ci  a  Timus-Qlsen  machine.  The  tip  masses 
were  remc  ed  and  the  lead  applied  to  the  twe  ends  simultaneously  by  an 
attachment  t j  the  moving  head  ’f  the  machine.  Moments  were  deduced 
from  strain  measurements  at  gauges  1  and  4  and  curvature  from  gauges 
2  and  3,  as  xr.  the  dynamic  tests-  Moments  were  also  calculated  from 
the  tccal  lead  applied  by  the  machine  These  were  not  in  close  agree¬ 
ment  with  the  mcments  measured  at  stations  i  and  4.  The  difference 
was  assumed  tc  be  due  tc  moment  restraint  provided  by  the  "hinged" 
supp.rts 

Material 

The  material  was  a  ’020  steel,  cold-rolled  and  annealed,  after 
machining,  at  800  C  icr  one  hour 
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Test  Results 

Measurements  from  gauges  1  and  4  showed  that  the  arrangement  and 
response  were  symmetrical*  It  was  assumed  that  the  moment  was  not  a 
function  of  position,  between  the  supports.  In  Figure  3  the  moment 
versus  curvature  from  two  static  tests  is  plotted  as  obtained  at  gauges 
2  and  3.  In  Figure  4  a  static  curve  is  compared  with  a  dynamic  curve 
obtained  from  a  10  inch  drop  of  the  shock  machine.  The  data  tends  to 
confirm  that  obtained  on  the  Naval  Research  Laboratory  shock  machine 
in  that  the  elastic  stiffness  El  is  •’•'out  the  same  under  dynamic  and 
static  loading.  Also,  shortly  che  outer  fibers  yield,  at  about 

120  inch-pounds,  the  dynamic  moment  becomes  less  thou  the  static 
moment  for  the  same  curvature.  As  the  loading  increases,  the  dynamic 
moment  again  exceeds  the  static  moment.  Thus,  the  "sag”  in  the 
dynamic  moment  curvature  curve  again  is  evident- 

In  Figure  5  the  same  dynamic  curve  shewn  in  Figure  4  is  shown  as 
curve  1  and  compared  with  that  from  another  specimen.  The  curve 
numbered  1  represented  the  first  test  of  a  beam-  The  curves  2,  3,  4 
and  5  were  successive  10  inch  drops  of  another  beam  that  had  previously 
been  dropped  several  times.  The  curves  show  the  effect  of  reloading. 
The  behavior  is  not  completely  understood,  but  this  trend  on  successive 
drops  toward  greater  curvature  for  the  same  moment  is  obvious. 

In  Figure  6,  >  he  result  of  increasing  the  drop  height  on 
successive  tests  is  shown-  Curve  1  is  the  result  of  a  10  inch  drop 
following  the  drop  which  produced  curve  5  of  Figure  5  Curves  2  and  3 
resulted  from  successive  12  and  14  inch  drops  Note  that  the  curvature 
scale  on  Figure  6  is  different  from  that  of  Figures  3,  4  and  5. 
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One  of  the  points  of  primary  interest  in  the  present  data  is  the 
time  delay  between  the  occurrence  of  maximum  moment  and  maximum 
curvature.  This  time  delay  can  be  directly  related  to  the  constant  C 
in  Equation  (1).  The  following  simple  mathematical  exercise  will  help 
show  this  relationship. 


Relationship  Between  C  and  the  Time  Delay  , 

The  purpose  here  is  to  assume  a  simple  form  of  moment  loading  and 
solve  Equation  (1)  for  the  resulting  curvature  in  order  to  relate  the 
constants  to  associated  phase  angle  or  time  delay  between  maximum 
moment  and  maximum  curvature. 

During  the  elastic  range  it  is  immaterial  as  to  the  time  and  rate 
of  loading.  In  the  plastic  range  the  loading  is  assumed  to  be 
sinusoidal.  The  following  theoretical  moment-time  input  is  used: 

M  =  S.  k  =  S.c  t  o£t<t  (2) 

1  1  o  o 

M  =  S.  k  =  M  t  =  to  (3) 

1  y  y 

Sc 

M  =  S.  k  +  - —  sin  w(t-to)  t>to  (43 

1  y  03 

Here  k^  is  the  curvature  at  the  knee  of  a  typical  moment-curvature  plot, 
and  is  attained  when  nearly  the  entire  cross-section  begins  to  undergo 
plastic  deformation.  Also  -  El  is  the  elastic  bending  stiffness. 

In  terras  of  symbols  used  in  reference  0),  Equation  (1)  could  be 


written 


K  +  R(M-M  ) 

S  L 


(5) 


Here 


R  SL  =  C. 


During  elastic  loading  M  =  M  ,  so 


k  =  c  t 
o 


o£t£to 


(6) 


and  k  =  k 


t  =  to 


(7) 


For  tito,  we  assume  a  bi-linear  static  moment  curvature  relation¬ 


ship: 


Mst  ■  M,  +  S2  <k'V 


(8) 


Substituting  (8)  into  (5)  yields  the  following  linear  differential 
equation: 


k  +  RS0k  =  •—  +  RM  +  R[S„k  -M  ]  for  Uto 
2  S1  2  y  yJ 


(9) 


With  M  as  prescribed  in  (5),  the  resulting  solution  for  k  is 
-(t-to) 


RS„ 


k  =  k  +  d. e 

y  1 


-  d^  c-'-s  w(t-to)  +  d2  sin  w(t-to)  (10) 


Whe  re 


dl  = 


R  a  w(S1-S2) 
Sx  (w2-RS2) 
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d2  = 


(aco2  +  R2aS1S2) 
S1(o)2  -  RS2) 


and 


S.c 
1  o 

w 


Considering  only  the  sine  and  cosine  terms  in  (10)  the  phase  angle 
(j>  between  moment  as  prescribed  in  (4)  and  curvature  k  is 

dl 

tan  <}>  =  -j— 


tan  <j>  = 


R  WCS^) 
_  - 

w  +  RS^ 


(11) 


2  2 
“  ✓  >n  r 


For  the  present  data,  it  appears  that  u)  <<R  S^S2 


Then  tan  <{>  = 


W(S2_S1) 


(12) 


If,  in  addition,  S2<<S^ 


tan  <J>  =  - 


0) 

RS, 


(13) 


Letting  <j»  =  WT,  with  T  the  time  delay, 
and  tan  <}>  =  <}>,  then  from  (12) 


T  = 


VS1 

RS1S2 


(14) 
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Since  C  =  RS^,  from  (14) 


C  = 


VS1 

TS„ 


(15) 


The  value  of  I  can  be  scaled  directly  from  data,  such  as  that  of 

_2 

Figure  7.  If  x  =  -0,2  x  10  sec.  and  =  0.2  S^, 

C  s  2000. 

With  S1  =  5000,  R  s  0.4. 

Note  that  RS2  £  1600  and  the  exponential  term  in  Equation  (10) 
drops  out  quickly. 


Results  and  Conclusions 

The  purpose  of  the  work  was  to  test  the  same  rectangular  beam 
cross-section  using  a  different  apparatus  for  direct  measurement  of  the 
moment-curvature  relationship  on  first  loading: 

(1)  The  result  of  previous  tests  with  accelerometers  appears  to  be 
confirmed;  namely,  that  the  dynamic  moment  becomes  less  than  the 
static  moment  shortly  after  yielding  of  the  outer  beam  fibers 
occurs.  This  is  apparently  in  conflict  with  Equation  (1)  and  the 
common  assumption  that  plastic  flow  will  not  occur  unless  the 
dynamic  stress  is  greater  than  the  static  stress. 

(2)  An  equation  is  presented  for  estimating  the  value  of  C  in 
Equation  (1)  directly  from  the  measured  time  delay.  The  value 
obtained  for  the  present  data  for  commercially  annealed  steel  is 
in  the  range  of  that  proposed  by  Brown ^  and  Vogel ^  based  on 
more  tho  mgh  analyses. 


►  — ^  <*  <f*w*£^fcs<. , 


*^r»- 
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(3)  The  apparatus  has  several  advantages  over  a  cantilever  beam 
arrangement.  One  is  that  there  is  practically  no  shear  in  the 
specimen.  A  second  is  that  some  of  the  difficulties  of  deducing 
load  from  acceleration  measurements  are  eliminated.  A  third  is 
that  the  distance  between  the  point  where  moment  is  sensed  and 
curvature  measured  is  greatly  reduced.  Also  most  of  the  material 
between  the  two  points  remains  elastic,  so  there  is  little  like¬ 
lihood  of  the  time  delay  being  due  to  wave  -travel  time  between  the 
two  points. 

(4)  The  study  has  dealt  primarily  with  first  loading  of  the  speci¬ 
men.  However,  the  data  from  reloading  of  the  same  specimen, 
especially  Figure  5,  bears  further  study.  Brown  gave  examples  of 
similar  behavior.  The  fact  that  the  moment  is  always  less  on 
reloading  than  on  the  previous  loading  was  initially  unexpected. 

There  had  been  speculation  that  Equation  (1)  could  be  used  for 
th 

the  i  1  reloading  if  M  were  replaced  by  the  relationship 

obtained  from  the  previous  dynamic  loading.  The  present  data 
indicates  that  prediction  of  response  to  reloading  requires  con¬ 
siderably  more  study. 


Fig.  7  Strain  Records  vs.  Time  on  Oscilloscope 
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II.  Relationship  Between  Moment-Curvature  and 
Stress-Strain  Constitutive  Relationships 

W.  Vogel  and  V.  H.  Neubert 


Much  basic:  study  has  been  devoted  to  the  behavior  of  materials  in 
simple  tension  or  compression.  Tests  have  been  performed  which  have 
led  to  constitutive  relationships  including  rate  effects.  The  purpose 
here  is  to  attempt  to  develop  the  moment-curvature  relationship  for  a 
beam  based  on  an  assumed  stress-strain  relationship.  The  development 

follows  that  presented  first  by  Vogel'  .  We  assume  that  the  total 

•  • 

strain  rate  e  could  be  represented  as  the  sum  of  e  the  elastic  strain 

rate,  and  e",  the  plastic  strain  rate,  that  is 

e  =  e'  +  e".  (16) 


In  Figure  8  a  cantilever  beam  carrying  a  concentrated  load  is 
shown.  If  the  beam  is  long  and  slender,  it  is  assumed  that  shear 
effects  may  be  neglected  with  regard  to  effect  on  yielding. 

Vogelv  '  assumed  that  the  total  stress  o  in  the  plastic  region,  of 
area  could  be  represented  by 

0=0+  be"  (17) 

St 


The  moment  M 

P 


due  to  stress  in  the  plastic  region,  if  the  beam  is 
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symmetrical  with  respect  to  the  neutral  axis  is 


M  =  /.  a  w(y)ydy  (18) 

p  a2 

where  w(y)  is  the  width  of  the  cross-section. 

Using  (17)  in  (18) 

M  =  /  a  w(y)ydy  +  b  /  e"  w(y)ydy  (19) 

p  s  c 

The  first  integral  is  the  contribution  to  M  due  to  the  stresses  un 

•  •  • 

A2-  In  the  second  integral,  let  e"  »  £  -  s'  from  (16). 


M 


’  (MstV 


+  b 


/  (£-£')  w(y)ydy 

A2 


(20) 


•  « 

If  plane  sections  before  bending  remain  plane  after  Vending  £  ~  ky, 

* 

with  k  the  curvature  at  the  cross- section,  Also,,  we  assume  £'  = 
and  (20)  becomes 


b  (M) . 


M  -  (M  ).  +  b  k  I. - =r 

p  st  A2  A2  E 


(21) 


On  the  elastic  region  0  =  ^st  so  contribution  to  the  moment 

M  due  to  stresses  on  A„  would  be 
s  2 


“s  -  <Va, 


The  total  moment  is  M  =  M 

P 


H*  M  ,  or 
s  * 


(22) 


rej  Q 


jCSW? 
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b  M. 


or 


M  =  Mst  +  b  k  h2  E 

Ma 

*  =  wr  +  bir  (M  -  V 


(23) 


(24) 


As  the  bending  increases  A ^  A  as  A^  ■+■  o  where  A  is  the  total 

area  of  the  cross-section,  resulting  in 


*  -  W  +  bT  ®  -  V 


(25) 


Although  there  are  assumptions  in  the  derivation  of  Equation  (24)  which 
are  open  to  question,  the  form  has  some  logic.  Written  in  the  form 


(M  -  Mst) 


(26) 


it  is  clear  that  when  the  section  is  entirely  elastic,  I.  =  o  and 

z 

M  =  o  resulting  in  M  =  M  . 

Comparing  Equation  (25)  with  Equation  (1)  and  Equation  (5)  in 
Part  I,  it  is  clear  that  C  =  —  =  RE1.  If  this  xs  true,  the  results  of 
the  present  study  may  De  extended  directly  to  beams  of  other  cross- 
sections  without  further  tests.  An  investigation  is  in  progress  to 
check  this  result. 

Equation  (25)  may  be  derived  in  a  more  direct  way  if  one  uses  an 

(4) 

equation  presented  by  Malvern  as  follows: 


(27) 
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Since  this  applies  to  the  complete  loading  range,  we  need  not  distin- 

•  • 

guish  between  elastic  and  plastic  regions.  Letting  £  =  ky,  we  multiply 
(27)  by  y  and  integrate  over  the  cross-section,  arriving  directly  at 
the  relationship 

!k  =  |  +  1  (28) 
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III.  Theoretical  Response  of  an  Idealized 
Viscoplastic  Cantilever  Beam 

R.  Weiss  and  V.  H.  Neubert 


Introduction 

One  of  the  goals  of  the  present  work  is  to  arrive  at  equations 
suitable  for  design  analysis  of  structures  loaded  dynamically  into  the 
plastic  range.  For  many  existing  computer  programs  using  finite 
elements,  it  is  easy  to  include  a  new  element  once  its  behavior  is 
stated  in  the  form  of  a  transfer  matrix,  a  flexibility  matrix,  or  a 
stiffness  matrix.  With  this  in  mind,  the  transfer  matrix  and  the 
stiffness  matrix  were  presented  in  reference  (1)  for  a  viscoplastic 
beam  element  which  behaves  according  to  the  equation  discussed  in 
Part'.-;  I  and  II  of  this  report: 

k  «  -=y  +  R  (M  -  M  )  (29) 

El  St 

To  retain  a  linear  form  of  the  differential  equation  it  was  assumed 

that  M  could  be  represented  in  a  piecewise  linear  manner,  the 
s  t 

simplest  form  being 


M  =  S,  k  o^-k^k  (30) 

st  1  y 


M  =  M  +  S0  (k  -  k  )  k.ik 
st  y  2  y  y 
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which  is  shown  in  Figure  9.  A  brief  discussion  was  presented  in 
reference  (1)  of  the  form  of  solution  of  the  differential  equation  of 
a  cantilever  beam  under  dynamic  loading  assuming  that  Equation  (29) 
applies  for  the  entire  beam  length. 

To  gain  experience  with  the  analysis  based  on  the  matrix  equations 
presented  in  reference  (1),  which  are  valid  in  both  elastic  and  plastic 
ranges,  a  computer  program  has  been  written  for  the  exact  solution  of 
a  two-segment  cantilever  carrying  a  tip  mass  as  shown  in  Figure  10. 

The  need  for  two  segments  arises  from  the  two  ranges  of  curvature  in 
the  bi-linear  form  of  M  t  as  given  in  Equation  (30).  In  order  to 
account  for  the  icn  beam  segment  in  the  constitutive  relationship,  an 
additional  subscript  is  now  inserted,  so  that  Equation  (30)  became 


c£ks.lc 

y 


M=M  +  S.0  (k  -  k  )  kkk 

st  y  i2  y  y 


(31) 


The  first  subscript  is  the  beam  segment  number,  the  second  subscript 
the  loading  range  on  the  linearized  M  versus  k  range.  With  the  bi¬ 
linear  curve  used  here,  implies  the  elastic  range  and  the 

plastic  range. 

The  loading  to  the  beam  of  Figure  10  was  taken  as  a  base  acceler¬ 
ation  yo(t)  of  a  triangular  form  as  shown  in  Figure  11.  This  approxi¬ 
mates  the  pulse  obtained  from  the  Barry  Shock  Machine  at  the  Naval 
Research  Laboratory  as  discussed  in  reference  (!)• 


Fig.  9  Bi-linear  Static  Moment  vs.  k 


Fig.  10  Two-Segment  Viscoplastic  Beam 


sj? 


-**—*+**»  —  *•«»&•«* 
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Fig.  11  yo  vs.  t 


As  the  beam  is  initially  loaded  =  L  the  total  beam  length  and 
l>2  =  0.  As  the  loading  increases  the  moment  at  the  base  would  first 
exceed  and  a  plastic  "hinge"  begin  to  develop.  If  at  this  time  a 
plastic  hinge  is  inserted  having  a  finite  length  I^,  the  moment  at  the 
elastic-plastic  interface  between  beam  segments  1  and  2  would  be 
slightly  less  than  M^,  but  would  continue  to  increase  as  the  dynamic 
loading  increases.  As  soon  as  the  moment  at  the  interface  equals  M  , 
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the  length  would  be  increased.  The  rate  at  which  the  interface 
moves  depends  on  the  material  and  the  loading  rate.  To  gain  experience 
with  the  effect  of  hinge  length  on  beam  response,  hinges  of  various 
lengths  were  used  and  their  lengths  keep  constant  throughout  the 
analysis.  This  is  equivalent  to  analyzing  a  beam  made  of  two  materials 
where  the  material  in  segment  1  has  such  a  high  yield  that  it  always 
remains  elastic.  In  later  analysis  this  limitation  is  removed  so  that 
L2  =  L2(t). 


Equations  of  Motion 

The  equation  of  motion  of  the  elastic-plastic  cantilever  beam  in 
terms  of  displacement  y^  of  the  tip  mass  relative  to  the  base  is 


•  •  •  •  » 


•  •  •  •  • 


A1  yl  +  A2  yl  +  A3  yl  +  A4  yl  =  B1  yo  +  B2  yo  +  B3^"My  +  S22 


(32) 


where 


Ai- 1 

D  A2  =  q  R  S21S22  +  31.2  L.  ,  SUS21  +  R  SnS21  +  q  R 


3  SUS n 

D  A,  =  - 

3  m. 


DA,  = 
4 


3R  S11S21S22 

m. 
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D  B. 


3L2  (4  +  L2/2)  R  SuS21 

m. 


B3  =  B3  (-  My  +  S22V 


D  ■  4  S21  +  3ll  h  S11  +  3L1  L2  S11  +  4  S11 


For  the  prescribed  pulse,  the  solution  for  y^  can  be  expressed 
exactly  during  the  three  stages  of  loading  in  Figure  11,  the  rise  era, 
the  decay  era,  and  the  free  vibration  era.  The  solutions  written 
below  are  in  terms  of  generalized  initial  conditions  occurring  at  yield 
time  t  at  the  onset  of  development  of  the  plastic  hinge. 

For  the  rise  era,  o^t^t^,  the  initial  conditions  are 


*1  (t  ’  V  '  *lr 

*1  ‘  V  ‘  hr  (33) 

*1  (t  ’  V  ‘  hr 


Hie  solution  is 


yi(t)  -  cx  "  V  +  ^  lc2  cos  r3(t  "  V1 


(y  )  max 

+  c3  sln  r3(t  ‘  V  +  TT~ 
3  4  r 


B1  +  V 


A3B2 


* 

B3 


(34) 
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1 

C1  2 
r 

s 


r2  +  r3  Glr  '  2r2  G2r  +  ylr| 


C2  "  r2  Vl(rl  ‘  2r2J  Glr  +  2r2  G2r  "  ylr 


-1 
-  y, ..  / 


C3  2  ^  rl 

r„r 
3  s 


2  2  \  „ 
r2  -  r3  -  Glr  + 


2  2  .  2\„ 

rl  ”  r2  +  r3  J  ^2r 


+  (r2  -  V  ylr 


2  ,  x  2  2 

rs  =  (rl  "  r2}  +  r3 


Glr  ylr 


G2r  "  ylr 


(y  )  max 
Jo 

A.  t 

4  r 


(yQ)  max 

A.  t 
4  r 


A-s 

B,  +  B.  t  -  -4- — 

1  2  y  A/, 


For  the  decay  era,  the  initial  conditions  are  used  in  the  form  of 
(33),  except  the  subscript  is  d  instead  of  r.  The  solution  is: 


rl(t  "  tv)  r2(t  “  tv) 

y^t)  =  D^e  y  +  e  y  [D2  cos  r,j(t  -  t  ) 


(y  )  max 

+  °3  Sin  r3(t  ~  V1  -  A~(td  -  tr) 


B1  "  B2(td  “ 


A3B2 


(35) 
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D1  =  ~T 

r 

s 


2  2 
r2  +  r3 


Gld  "  2r2  G2d  +  yid l 


D2  =  ^2  Vl(rI  "  2r2^  Gld  +  2r2  G2d  "  yld  | 

s 


D  -  1  ■( 

f 

2 

2  \ 

2 

2  u.  2 

°3  2  i 

r3  rs 

Lri) 

r2  ' 

'  3  12 |  Gld  + 

1 

rl  ' 

'  r2  +  r3 

2d 


+  (r2  -  ri>  yu 


(y  )  max 

Gld  “  yld  +  A. (t ,  -  t  ) 
4  d  r 


A„B„ 

B.  -  B,  (t  -  t  )  - 
1  2  d  y  A, 


A. 

4 


.  (y&)  max  B2 

G2d  =  yld  +  A. (t  '  -  t  ) 
4  d  r 


After  the  pulse  is  terminated,  for  t^t^,  the  initial  conditions 
are  used  as  in  (33)  with  an  a  subscript  instead  of  an  r.  The  solutio*n 


is : 


y^(t)  =  E1  e 


x  r2(t  ”  ty) 


1'  y  +  e 


[E2  cos  r3(t  -  ty) 


+  E_  sin  r,(t  -  t  ) ]  + 

3  3  y 


(36) 
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To  more  effectively  study  the  effects  of  key  parameters  such  as 
the  material  constant  R  and  plastic  hinge  length  L a  computer  program 
was  written  which  calculates  tip  displacement,  velocity  and  acceler¬ 
ation,  base  moment  and  curvature,  and  the  time  rates  of  change  of  both 
base  moment  and  base  curvature  at  various  intervals  of  time  until  the 
dynamic  moment  at  the  base  equals  the  yield  value.  When  the  beam  has 
yielded  at  the  base,  a  plastic  hinge  is  inserted  at  time  t^,  and  again 
the  response  variables  calculated,  including  the  moment  at  the  elastic- 
plastic  interface.  Realistic  values  for  the  base  acceleration  were 
obtained  from  reference  (1),  Figure  11,  page  110.  The  values  used  are: 

(y  )  max  =  100  g  m.  =  0.012  //sec^/in 

o  i 

2 

t  =  0.004  sec.  Sn ,  =  S01  =  5000  //in 

r  11  21 

2 

t^  =  0.006  sec.  S2J,  =  500  it  in 


L  =  =  5  inches 


1 


31 

Values  of  yield  moment  (247  ini?)  and  yield  curvature  (.05  in  *)  were 
gotten  from  reference  (3). 

In  Figures  12-16,  the  tip  displacement  versus  time  is  shown  for 
various  values  of  R  and  L ^  It  was  hoped  that  some  of  the  response 
variables  might  be  found  to  be  relatively  insensitive  to  hinge  length 
and  that  this  would  simplify  design  analysis.  To  date,  it  has  been 
found  that  hinge  length  is  quite  important.  On  each  figure  one  example 
curve  is  also  shown  for  R  =  1  in  which  the  hinge  was  inserted  at  t  =  0. 
The  effect  on  prediction  of  maximum  displacement  is  an  error  of  20%. 

In  Figures  17-20,  the  moment  M2  at  the  elastic-plastic  interface  is 
plotted  versus  time.  As  soon  as  this  moment  reaches  the  yield  moment 
the  hinge  would  be  extended  if  both  beam  segments  were  of  the  same 
value.  The  main  use  for  the  curves  at  present  is  that  they  give  an 
analyst  an  idea  of  the  rate  of  propogation  of  a  hinge  for  various 
values  of  L9  and  R. 

Results  Using  a  Variable 

The  next  step  was  to  use  a  variable  viscoplastic  hinge  length  I^. 
For  the  present  examples,  1^  was  incremented  in  0.05  inch  intervals. 

The  process  was  as  follows.  As  in  the  previous  examples,  the  moment 
at  the  support  was  monitored.  As  soon  as  it  reached  the  yield  moment 
M  ,  a  viscoplastic  hinge  0.05  inches  long  was  inserted.  The  moment  M„ 

y  *• 

at  the  elastic-viscoplastic  interface  would  at  that  instant  be  slightly 
less  than  M  .  The  timewise  variation  in  M„  was  then  monitored  and,  as 
soon  as  M2  =  M^,  the  hinge  was  lengthened  by  .05  inches.  This  process 
was  continued  to  allow  the  hinge  to  develop  as  required  by  conditions 
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at  the  interface.  The  effect  of  the  incrementing  procedure  and  its 
relation  to  convergence  is  being  studied  further. 

In  Figure  21,  the  effect  on  the  tip  displacement  y^(t)  of  having 

a  variable  hinge  len6 _  is  indicated  for  R  =  0.1.  The  solution  for 

1*2  a  variable  is  compared  with  those  previously  discussed  for  =  0.05, 
0.25,  and  0.50  inches.  It  can  be  seen  that  the  variable  hinge  solution 
is  very  close  to  that  for  =  0.50  inches.  In  Figure  22  a  similar  set 
of  curves  is  shown  for  R  =  1.  Here  the  more  exact  variable  hinge 
solution  is  closest  to  that  for  =  0.25  inches. 

In  Figure  23  the  variation  of  with  time  is  plotted  for  R  =  0.1 
and  1.  For  R  =  0.1,  the  hinge  reaches  its  final  length  of  0.5  inches 
in  a  relatively  short  time.  For  R  =  1,  the  hinge  spreads  less  rapidly 
and  reaches  a  total  length  of  0.3  inches.  It  is  interesting  to  note 
that  these  total  hinge  lengths  are  approximately  the  same  as  the 
constant  lengths  whose  solutions  for  y^  on  Figures  21  and  22  most 
closely  agreed  with  those  for  variable  hinge  length  as  discussed  in 
the  previous  paragraph.  This  may  have  some  practical  utility,  since 
it  is  easier  to  solve  for  response  .g  a  constant  hinge  length. 
However,  at  present  we  have  no  way  to  predict  the  hinge  length  without 
carrying  out  the  more  exact  solution. 

The  variation  of  moment  M2  with  time  at  the  moving  interface  is 
shown  in  Figure  24  for  R  =  0.1  and  1,  For  R  =  0.1,  the  hinge  develops 
to  its  full  length  rapidly  and  then  the  interface  moment  begins  to 
decrease.  For  R  =  1,  a  cyclic  decrease  and  subsequent  increase  in 
moment  each  time  the  length  of  the  hinge  is  incremented  is  much  more 
pronounced,  indicating  that  an  increment  length  smaller  than  0.05 


inches  should  be  used. 
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Summary  and  Conclusions 

The  purpose  of  this  part  of  the  report  has  been  ko  present  some 
results  from  a  computer  program  for  the  exact  analysis  of  the  visco¬ 
plastic  beam  under  shock.  The  displacements  calculated  using  hinges 
of  constant  length  and  variable  length  L ^  should  be  upper  bounds, 
because  elastic  unloading  and  shortening  of  the  hinge  have  not  been 
accounted  for  in  the  curves  presented.  The  computer  analysis  can 
easily  take  these  effects  into  account.  One  of  the  final  goals  of  this 
theoretical  study  is  to  investigate  the  effects  of  yielding  of  supports 
and  foundations  on  design  shock  spectra. 

In  the  overall  study  of  development  of  yielding  and  failure 
criteria,  the  importance  of  experimental  data  can  not  be  over¬ 
emphasized.  The  constitutive  relationship  postulated  in  Part  12 
accounting  for  dependence  on  cross-sectional  shape  of  a  beam  needs  to 
be  checked  experimentally,  and  modified,  if  necessary.  Different 
cross-sectional  shapes  of  various  practical  materials  need  to  b<. 
investigated.  In  general,  the  ability  to  calculate  using  the  digital 
computer  presently  exceeds  the  amount  of  experimental  data  available 
for  guidance  of  analysis. 


Fig.  12  y.  vs.  c,  Various  R,  L-  =  0.05  inches 


vs.  t,  Various  R,  L_  =  0.25  inches 


Fig.  14  y.  vs.  t,  Various  R,  L_  =  0.50  inches 


R  H  .  1 


Fig.  17  vs.  t,  Various  R,  L2  »  0.05  inches 


Yield  Moment 


Fig.  23  L,  vs.  t,  R  *  .1  and  1 
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